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Abstract 

We use the techniques of Bartnik fSI to show that the space of solutions to the Einstein- 
Yang-Mills constraint equations on an asymptotically flat manifold with one end and zero 
boundary components, has a Hilbert manifold structure; the Einstein-Maxwell system can be 
considered as a special case. This is equivalent to the property of linearisation stability, which 
was studied in depth throughout the 70s [T] [2l [9l [TT] [T3] [THl [1^ . 

This framework allows us to prove a conjecture of Sudarsky and Wald !22 , that is, the 
validity of the first law of black hole thermodynamics is a suitable condition for stationarity. 
Since we work with a single end and no boundary conditions, this is equivalent to critical points 
of the ADM mass subject to variations fixing the Yang-Mills charge corresponding exactly to 
stationary solutions. The natural extension to this work is to prove the second conjecture from 
[22], which is the case where an interior boundary is present; this will be addressed in future 
work. 

1 Introduction 

A solution to the full Einstein- Yang-Mills equations is given by a Lorentzian metric and a g— valued 
one- form on a 4-dimensional manifold M, where q is the Lie algebra of some compact Lie group 
G. It is well known, that given a sufficiently regular solution, (g, A, tt, e) of the constraint equations 
(|3.ip - (|3.3p on a 3-manifold, A4, we can find a full solution to the Einstein- Yang-Mills equations 
with an embedded hypersurface on which the initial data is induced. By this, we mean g is the 
induced metric; tt = [K — tr{K))y^, where K is the second fundamental form; A is the projection 
of the Yang-Mills connection onto the hypersurface; and e is four times the negative of the induced 
Yang-Mills electric field density, E, as viewed by a Gaussian normal set of observers. We use tt and e 
instead of K and E as these quantities are the canonical momenta for the Hamiltonian formulation 
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used. 



Geometrically, the Yang-Mills fields can be thought of as coming from some principal G-bundle, 
^P; for simplicity, we assume this bundle is trivial, "^P = '*A^ x G. The interpretation of the Yang- 
Mills potential, '^A is as the puUback of a connection one-form uj on via a global section. Given 
a global section l : — x G, we define '*A = l*uj and = 6*(ri), where H. = duj + uj A uj is 
the curvature form of uj. The Yang-Mills potential is a g- valued one- form on '^A4 and "^F is a 
0- valued two- form on , called the field strength tensor. The restriction of "^A and "^F to A4 can 
be viewed respectively, as the pullback via some section of a connection one-form and associated 
curvature on a bundle P = A4 x G, the restriction of to . 

The outline of this article is as follows: In section [2l we define a Hilbert manifold structure for 
the set of possible initial data, the phase space. In section [3l we use an implicit function theorem 
argument to prove that the set of solutions to the constraint equations, is a Hilbert submanifold of 
the phase space; we call this the constraint submanifold. In section |4] we define the energy, momen- 
tum and charge functionals on the phase space, and construct an appropriate Hamiltonian for this 
system, similar to that of Regge and Teitelboim [5^. This Hamiltonian ensures that Hamilton's 
equations give the correct evolution equations and on shell, it gives a value for the total energy of 
the system. In section [5l we use a Lagrange multiplier argument to prove that stationarity is not 
only a sufficient condition for the first law of black hole thermodynamics to hold, it is necessary. 
Evidence for this is given in [22] , however a rigorous proof requires the Hilbert manifold structure 
discussed in section [31 

The phase space considered is tuples (g,^,7r,e) with x x x local regularity, with 
appropriate decay conditions on g, tt for asymptotically flat spacetimes. The decay conditions on 
the fields A and e (discussed in section [2]) are more subtle; in addition to the requirement that 
the fields are asymptotically zero, we further require that A approaches a collection of Maxwell 
(photon) fields at a faster rate. For simplicity, we work on a 3-manifold Ai with one asymptotic 
end and no interior boundary, however it is clear that all results will still hold in the case of many 
asymptotic ends. We will consider separately the case when Ai has an interior boundary. 

The Hilbert manifold structure for the space of solutions is equivalent to the property of lineari- 
sation stability, which was studied by many authors throughout the 70s. Linearisation stability of 
Minkowski space was established by Choquet-Bruaht and Deser in 1973 [11] and in the same year, 
Fischer and Marsden proved linearisation stability for non-exceptionaQ data on a compact mani- 
fold pjj. This was subsequently extended to the Einstein-Maxwell |lj and Einstein- Yang-Mills j[2] 
cases by Arms. The general asymptotically flat case wasn't established until 2005 when Bartnik [5] 
provided a Hilbert manifold structure for the phase space for the Einstein equations. We follow the 
techniques of Bartnik to generalise this to the Einstein- Yang-Mills case. Like Bartnik, we consider 
a class of initial data too rough to guarantee that a solution to the constraints corresponds to a full 
solution, however linearisation stability can be obtained by noting that the analysis presented in 
section [3| remains valid if the phase space is required to have enough regularity for known existence 
and uniqueness theorems to be applied [T71 [5T] . 

^Moncrief later proved that the exceptional data considered by Fischer and Marsden corresponds exactly to 
solutions exhibiting Killing fields [181 119] 
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2 Notation and Preliminary Definitions 



Let be a paracompact, connected, oriented and non-compact 3-manifold without boundary, and 
suppose there exists a compact Mq C M and a diffeomorphism (p : A4 \ A4q ^ \ Bi{0) , where 
Bi{0) is the closed unit ball. Let g be the Lie algebra of some compact Lie group G, and recall 
that any such Lie algebra must be the direct sum of abelian and semi-simple Lie algebras. We can 
then define an adjoint invariant positive definite inner product, 7 on g, given by the negative of 
the Killing form on the semi-simple factor; we may use the regular Euclidean inner product on the 
abelian factor. 

Throughout this article, we use four different sets of indices on different objects as outlined below. 





Latin lower case, mid-alphabet ... 




Greek lower case, mid- alphabet 





Latin lower case, early alphabet a, 6... 




Greek lower case, early alphabet a,f3... 



Fix a smooth background metric g on A4, such that g = 0*(o-) on A4 \ Mq, where a is the 
Euclidian metric on M.^ . In terms of this background metric, we define the weighted Lebesgue (Lg) 
and Sobolev {W^''^) spaces as the completion of C^{M) with respect to the norms 

ikiu, = | (/A.i^r''"''"'^^")'^', P<oo (2.1) 

[ esssup_;^(r \u\), p = oo 

fc 

ll"IU,P,^ = El|V^'^^IU-„ (2.2) 

j=o 

where r > 1 is a smooth function on A4 and r{x) — \4>{x)\ on \ A4q is the regular Euclidian 
distance from the origin in M.^. Objects labeled with an 'o' refer to objects associated with g, such 
as the connection V and volume form dfj,^. Spaces of sections of bundles are defined on Ai in 
the regular way with respect to g and 7 where appropriate and we will omit specifying the bundle 
where there is no risk of confusion. Intuitively, the spaces defined above contain functions of local 
regularity or W^'^ , which behave as o{r^) as r — >■ 00, with derivatives decaying appropriately. 
Since we are working with a trivial bundle, we may choose the flat connection as a background and 
work in a gauge such that the background gauge covariant derivative is exactly V. 

The usual definition of Yang- Mills total charge [12] is given by 

4^Q" := lim / (2.3) 

where Sj^ :^ {x E Ai : r{x) = R} is the sphere of Euclidean radius, R, for large R. Unfortunately 
this is gauge dependent in general; in fact, it may be that the integral is finite in some gauge 
and infinite in another. A sufficient condition to ensure the charge is well defined, is to ensure 
[A^, F^"^] e (cf. 12 ). In the Hamiltonian formulation, the choice of Aq is still a gauge freedom, 
however we instead impose the condition that [ylj,£'"'] e for any possible choice of E in the 
phase space. 
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Specifically, we will enforce that the dominant part of A near infinity is valued in the centre of 
0. In the language of physics, this condition is the that the gauge field behaves like a collection of 
photon fields near infinity. In the language of physics, the non-commutativity of A corresponds to 
self-interactions which would cause the fields to drop off quickly; it is therefore not surprising that 
we require such a condition. 

In defining the following spaces, we will make use of the decomposition 9=3©?, where 3 is the 
centre of q. 

g : = W^!f/2 (^2) JC: = W^l^ (S^ ® A^) 

A : = T4^!f/2 iT*M 3)® 1^^3/2 (r*X (g) J) f : = W^^^^ [TM ® Q* ® K^) 

JV : = L^i/2(A" xTMx g(E)A°) Af* : = L^5/2(A^ x T*M ® A^ x g* (g) A^), 

where A*^ are fc-forms on A4 and 5*2 and S'^ are symmetric covariant and contravariant 2-tensors 
on A4 respectively. The spaces Af and Af* can be interpreted as spaces of vector fields and 
covector densities, respectively on ^P. The direct sum in the definition of A is understood as 
the internal sum in W_y^{T*A4 (E) g). For an arbitrary A ^ A^ + A, e A, we will write 
W^Wa ■— 11^3112,2,-1/2 + ll^{||2,2,-3/2- Note, if = u(l) = 3, then this includes the regular de- 
cay conditions {E,B — 0{r^'^)) for the Einstein-Maxwell equations. 

Define the spaces 

G'^ = {9 -9-9(^5,9 >0} 

g+ = {,g e g+ : Ag < g < A^^^} , < A < 1. 

From the weighted version of Morrey's inequality ()2.10p . we can deduce that both g E G and A E A 
are Holder continuous with exponent ^. In particular, the inequalities in the definitions of Q'^ and 
Gx are understood in the pointwise sense. 

The phase space we will consider is 

T -.= 0+ X Ax JCx £, (2.4) 



which is independent of g (and (f>) (cf. [5]). 

We now quote directly, the weighted Sobolev-type inequalities from [4]. 
Theorem 2.1. The following inequalities hold: 
i-) -(f 1 < P < 9 < 00, 62 < ^1 O'f^d u e L\^, then 

IklU <c|kll,,^, (2.5) 

and thus L^^ C . 
a.) (Holder) If u E Lg^, v G Lg,^ and d = Si + 62, ^ < q,r < 00, then 

h^Wp.S < Mks, (2-6) 

where 1/p = 1/q + 1/r. 
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Hi.) (Interpolation) For any e > 0, there is a C(e) such that, for all u e Wj'^ 

IHIi,p,5<e|l«lkp,5 + C(e)N|,,,, (2.7) 

for 1 < p < q < oo. 
iv.) (Sobolev) IfueW^'^, then 

\Mnpnn-kp),S<^Mk,,,6 (2-8) 

for q satisfying p < q < np/{n — kp). 
If kp> n then 

\Moo,S < c||M||fe,p,5 (2.9) 
V.) (Morrey's) If u G W^'^ and 0<a<k — n/p<l, then 

< 4'4k,p,s, (2.10) 

where the weighted Holder norm is given by 

||u|Uo,a := sup (r (x) sup — j ^ 

^ xeM^ 4\x-y\<r{x) F - 2/1 



sup 
xeM 



{r\x)\u{x)\j 



vi.) (Poincare) If 6 < and 1 <p < oo, for any u G W^'^ we have 

\\u\\p,s < c||Vu||p,5_i, (2.11) 

where n = 3 is the dimension of A4 . 

3 The Constraint Submanifold 

The constraint equations - defining the constraint map, for sufficiently smooth data - are given 

by 

^o{g,A,7r,s) = (^(Tr^)' - tt'^tt,^ - 2{E^,El + BtBt))g-'/^ + 

= ili^f - ^''^ij - ile'ael + 2BlBl))g-"^ + R^g = Too (3-1) 
$,(5, A, TT, 6) = 2V^n,, - ei{V,A] - VjA^ + CtA^^i) = Toi (3.2) 

{g, A, n, e) = -djsi - C^.A^jsi = (3.3) 

where (r^oiio) is some prescribed source. The quantity := ^e^'^{djAg^f. — d/^A^j + C^i,AjAl) is 
the Yang-Mills magnetic field, as viewed by a Gaussian normal set of observers; C^^ ~ C^bc ^-re the 
structure constants of q and e*-'*^ (resp. ej,^.) is the completely antisymmetric tensor density with 
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weight 1 (resp. —1). Also note, since e is a vector density, we have V ■ e — d ■ e — V ■ e. 
It should be noted that the quantity e = ~AE differs from the canonical momentum of Arms [2] by 
a factor of 4, and is the negative of that used by Sudarsky and Wald [22]; we will alternate between 
using E and e, wherever it is convenient. As usual, we have used natural units where c — G — \ and 
we have set the coupling constant to -\/47r to agree with regular electromagnetic theory; specifically, 
our constraints are derived from the action 



5 = / [R-W]"). 

J M 



(3.4) 



Note, this differs from the regular action by a factor of IGtt. If we were to use a different coupling 
constant, a factor would be present in the term, however this makes no difference to the argu- 
ments presented. 



Since we are following arguments from [5], it will be useful to define the pure gravitational 
constraint map 



*(5,7r) := 



$0(3, A, Ti, e) + 2{e':e1, + B^,Bt)g-^^^ 









*»(.9,7r) 



(3.5) 



Throughout, we use c or C to denote some constant depending on {Ai^g) and other fixed param- 
eters, which may vary from line to line. Where appropriate, we will make explicit the parameters 
on which these constants depend. 

We first show $ : — > M* is a smooth map of Hilbert manifolds. 

Proposition 3.1. Suppose (5,^, 7r,e) e Q'l xAxJCxScT for some fixed A > 0, then there 
exists a constant c = c(A) such that 

ll'^'o(ff; ^)ll2,-5/2 - ^(^ + 11.9 ^ .9112,2,-1/2 + ll^ll 1,2,-3/2 + ll^ll 1,2,-3/2 + 

(3.6) 

||$,,(5, A,^,e)|l2 -5/2 < c(||V7r||2,_5/2 + l|Vff||i,2,-3/2 Iklli 2 -3/2 + ll^lll 2 -3/2 (1 + WMI 2 -I/2)) 

(3.7) 

||$J5, A^,e)ll2,-5/2 <c\\e\\, (1 + PIU) (3.8) 
Proof. From [5] (Prop. 3.1) we have the bounds. 



1^0(5,^)11 



2,-5/ 



ll*»(5,^)ll2,-5/2 < c(||V7r||2^„5/2 

thus we need only to bound the Yang-Mills terms. 
Applying (|2.6p and (|2.8p . we have the inequality 



2 < C(l + II5 - 5ll2,2,-l/2 + ll'^lll,2,-3/2) 

^5lll,2,-3/2 lklll,2,-3/2 



(3.9) 
(3.10) 



h'^h.S < \\u\\l.5/2 < c||u||i,2,5/2- 



(3.11) 
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Using ([XTT]) . (HH) and ([2^, we have 

l|i3'll2,-5/2 < c(||Vv4||i,2.-3/2 + II [A ^] lk2,-3/2)' 

< c(Pll2,2,-l/2 + P{||l.2,-3/2)^ 

< c(P||2,2,-l/2 + P?ll2.-3/2 + l|V(At)A(||2,_5/2)' 

< C(l + ||AJ|2,2,-l/2 + Pfll?,2,-3/2 + l|VA(||4,_5/4||^||4^_5/4f 

< C(l + ||AJ|2,2^_i/2 + Pt||?,2,-3/2 + l|Vv4(||i,2,-5/2Ptlll,2,-3/2)' 

<c{l + \\A\\t,} 

The term is clearly taken care of by p. lip ; combining these bounds with the definition of , 
we have estabhshed p.6p . 

Similarly, we have 

||$,(5, A,^,e) - *,(5,7r)||2^_5/2 < c(||£V^||2,-5/2 + l|£4ll2,-5/2) (3.12) 

< c(||£||4,-5/4l|VA||4,-5/4 + l|e||4,-3/2ll All8,-l/2) (3.13) 

<c||e||l,2,-3/2(l + Pll2,2.-l/2), (3.14) 

which establishes p.7p . 
Finally we have 

||$,(g,A,^,e)||2,_5/2 <c(||£||i^2,-3/2 + Ptll4,-5/4l|e||4,-5/4) (3.15) 

< c(||£||l.2,-3/2 + Ptlll.2,-3/2l|£||l,2 -3/2) (3.16) 

giving p.Sp . and thus completing the proof. 

□ 

Corollary 3.2. $ : J" ^ TV* is smooth. 

Proof. It can be seen from Proposition 13. 1[ that $ : — )■ A/"* is locally bounded. We note, R can 
be expressed as a polynomial function in g, g~^, Vt; and vg, and therefore the constraint map 
can be considered a polynomial function in 12 variables, 

Hg, Va, 1/V5, V5, tt, Vtt, e, Ve, A, VA) = $(5, A, tt, e). (3.17) 

For positive definite matrices, the maps g 1— >■ Vg, g 1— >■ V^i?, A 1— >■ V>1, g 1— >■ etc. are smooth. 
Further, locally bounded polynomial functions are smooth (in the sense of Frechet differentiability) 
(cf. [IB], chapter 26), it follows that $ is a smooth map of Hilbert manifolds. 

□ 

The linearisation of $ at a point G = {g, A, TT,e) E J- 

D^G -.GxAxJCxS^M* 
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is given by 

D%aih, b,p, f) = (tt^^^- - 27:'\i - 2{ElE''= + BlB''=))h,^g-"' 



^2 4^ 

,1. 

-2' 



+ [pl-K] - 2^^-'p,,)g-i/2 - \raelg-^l^ (3.18) 

+ 2V,p^' ~ (V, A," - V, A" + Cfc"cA^A,=) (3.19) 
D^^ciK b,p, /) = - (C„V'c6' + V,/i + C^fc/^'A-). (3.20) 
See [M] (and references therein) for computations. 

The adjoint is simply computed by integration by parts and throwing out the boundary terms: 



D'S>;aiN,X,V) ^NU^kT^'^ - 27t'\1 - 2{ElE''' + BlB"') 



1/2 



- \{4f + {E^El + BlBl)^ g^^)f 
N{]^Rg'' - R'') + WN - V'^ VfcTvj ^ + Cx^"' (3.21) 



+ 



- X^(V,ei + C:,A]ei) + X^C^^A.^^ - C'^aslV' (3.22) 
i?a>;G(7V,X,F) =iV(g,^-7r,^ - 271,^ )5-1/2 _ (3 23) 

i?<i>:G(7V, X,V)^~ \Ne1g~"^ + ^A'] - V,A? + CIA\A'^) + a.V'^ + CIA\v' 

= NEtg'"' - e,,fcX^S"'- + + (3.24) 

The tuple {N, X,V) € N corresponds to a scalar function, vector field and g- valued function on M. 
respectively. We will omit reference to the base point G when there is no risk of confusion. 



It should be noted that the map given by 

T{N,X) := 



D^;{N,X,V) 
D^l {N,X, V) 



(3.25) 



is of the exact same form considered in [5], so we quote the following two theorems. 



Theorem 3.3. //(A, /a) G x W_y^ and {N,X) e L_i/2 is a weak solution ofT{N,X) 

(/11/3); then {N,X) G W'^f^j^ is a strong solution. 



Theorem 3.4. The operator T has trivial kernel in iTl 



1/2- 
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From this we will prove also has trivial kernel, but first wc will need to prove an estimate. 
Lemma 3.5. If ^ = iN,X,V) £ W^'^^^ satisfies D^l{£_) = f^e then 

11^112,2,-1/2 < c(||(iV,X)||2,2,o + ||/4|ll.2,-3/2 + ll^ll2,o), (3.26) 

where C depends on (g, A, tt, e). 
Proof. From (|3.24p we have 

= D^liC) - NEtg'^'^ + e.^kX'B"^ - CtA^V" (3.27) 

For brevity, let C = {N,X) £ W^'f^^ and ^ = {E,B) £ W^'^^^. By differentiating (p:?7|) and 
applying the inequalities from Theorem 12. 11 we have 

l|VV||2,_5/2 < c(||/4||l,2-3/2 + 1 1 V (C)/? 1 1 2,-5/2 + II V(/3)C|| 2,-5/2 
+ |V(Aj)t/||2^_5/2 + ||V(mil2,-5/2) 

< c(| /4||l,2-3/2 + l|V(C)||4,-lll/3||4,-3/2 + 1 1 V (/?) 1 1 2,-5/2 1 1 Cll oo,0 
+ I V^||4,-5/2ll^ll4,0 + l|Vl^||2,-l||A(||^._3/2) 

< c(| /4lll,2-3/2 + l|V(C)||l,2,-lll/5||l,2,-3/2 + 1 1 1 1 1,2,-3/2 II Cll 2,2,0 
+ I VA.||i,2,-5/2ll^lll,2,0 + l|VV||2,-l||A(||2,2,-3/2) 

< c(^| /4lll,2-3/2 + ll/?lll,2,-3/2llCll2,2,0 + 1 1 ^1 1 1,2,0 IIA 1 1 2,2,-3/2) 
<c(||/4||l,2-3/2+IICll2,2,0+ll't^lll,2,o), 

where the constant, c, depends on (/?, A^) in the last line. Applying the weighted Poincare inequality 
(PTT|) . we have 

l|l^ll2,2,-l/2 <c(||(Af,X)||2,2,0+ll/4lll,2,-3/2 + ll^^lll,2,o)- (3-28) 

Applying the interpolation inequality (j2.7p to the last term on the right hand side and choosing e 
small enough, gives us (|3.26p . 

□ 

Theorem 3.6. If £, £ £^1/2 is a weak solution of D<P* {£_) = {fi, f2, fs, h) , where {fiJsJi) £ 
L-3/2 X T^!3% X W'_'^/2' then ^ £ W^'^^^. 

Proof From Theorem [331 we have {N,X) £ W^Jy^. From it can be seen that on each 

VI GG M.,V (weakly) satisfies an equation of the form 

d^V- = at^^' + Pt (3.29) 

with coefficients a £ W'^''^, j3 £ W^''^ on ft. The argument here is well known; V can be approxi- 
mated by Vg, where e is a mollification parameter and on 51 we have 

IIVKII2 < c(||a||^||yj|2 + \\l3h) < c(||a||2,2ni2 + IWh)- (3.30) 
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Since is uniformly bounded in PF^'^, it follows V G W^'^ on il. 

By differentiating p.24p . it can be seen that on any fl CC A4, V weakly satisfies an equation of 
the form 

4- = ald^ V' + (3%cV' + e^^j , (3.31) 

2 2 1 2 2 1 2 

where here we have a G W ' , (5 G W ' , G L . Since we now have V G W ' , the same argument 
will give us — > V" e PF^'^. All that remains to show, is that V and its weak derivatives satisfy 
the correct asymptotic conditions. 

For any cutoff function xb. with x_r = 1 on i?fl(0) and zero outside B2fi{0), £, = {N,X,XrV) 
satisfies the conditions for lemma [3?5] and thus we have 

IIXfl^ll2,2,-l/2 < t^(ll(A^,^)ll2,2,-l/2 + ll/4lll,2,-3/2 + ll^llo,2,o)- (3-32) 

Once more we have a uniform bound, thus it follows Xr^ V in 

□ 

Next we demonstrate that D$* has trivial kernel. This amounts to saying that there are no 
symmetries of the data, asymptotic to zero at infinity - this will be discussed in more detail in 
section 21 

Proposition 3.7. If £, e satisfies D^*{£_) = on M, then f = 0. 

Proof. From theorems 13.41 and 13.61 respectively. {N^X) = and V G W'^^^. We have djV'^ = 

ClV'Al from (IX^ . so we can repeat the arguments with V^, as we easily obtain the uniform 
bound 

ll^j3,2,-l/2<c||5V,|l2,-7/2 

< c(||V%||2._7/2nioc.O + llVAf |U^_5/2|| W||4,_i + || A( |U^_3/2 || V V|l 2, ^2) 

< c||A(||2, 2 -3/211^^112, 2,0- 

As above, we now have V G W^f^j^; from the weighted version of Morrey's inequality (|2.10p . we 

1 1/2 

have V G C . That is V strongly satisfies the equation 

\d,{r'V^) = - Ctyj^V, - 0, (3.33) 

and since V is asymptotic to zero (and is connected), V = Q. 

□ 

We are now able to apply the implicit function theorem to prove the level sets of $ are smooth 
submanifolds of J^. 

Theorem 3.8. For any {s, Si,aa) G Af* , the set 

C(s, 5, a) {(<?, A,7T,e)GT: $(3, A, n, e) = (s, S, a)} (3.34) 
is a Hubert suhmanifold of J- . 
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Proof. We simply must establish that is surjective and splits its domain into the direct sum 
of the Ker{D^) and a complementary subspace, then the result follows from the implicit function 
theorem. Since is bounded, the kernel is closed and hence splits. The codomain splits as Af* — 
Ran(L'"I>) © coker(Z?$), but from proposition l3.71 coker(L'$) is trivial. To establish surjectivity, we 
simply must show that has closed range. 
Consider variations of the form 

= -^grjV ^1 = (3.35) 

^ 1 (v^r^- + v^r^" - Vfer'=5'^")V5 /f = -d^rVa- (3.36) 

With {h,b,pj) of this form, define F{y) = F{y,Y,'il;) = £'$(g,A,,r,£)(^, 0,p, /). Explicitly, 
F{y,Y,^) = 

Clearly we have F : W^'^^^ ~^ ■ This new operator is clearly bounded and the adjoint map has 
similar structure. We have the following scale-broken estimate from [3] 

ll"ll2,2,-i/2 < c(||Au||2 -5/2 + llujb.o), (3.37) 
from which we can establish an elliptic estimate for F. 

l|A3^||2,-5/2 < c(||F(J^)||2,-5/2 + ll*0y|l2.-5/2 + I KV J^l I 2, -5/2 + lkfy||2,_5/2 + ) 1 1 2, _5/2 

+ ||V(V')VA||2,_5/2 + l|V(^)A?||2,_5/2 + l|V(7r)y||2^_5/2 + \\<7{mth^5/2 
+ \\E^y\k-5/2 + l|S'y||2,-5/2 + II^Vl^ll2,-5/2) (3.38) 

Where f := f *j — r*^ = \g^^ ^ i9ji + ^ j9ii — '^i9ij) is the connection difference tensor and is 
clearly W^'^^r^. It's easy to check that Ric is of the form Ric ^ {Ric + Vf + f ^), so it follows 
Ric G 

For brevity, define Ui := ($o, -R«c, Vtt, 7rf , B^) e £^5/2 and U2 := (tt, V(A), A(, A( , i;) e 

1 2 

^-3/2- With this notation we have 

113^112,2,-1/2 < c(||F(3^)||2,-5/2 + \\Uiyh^-5/2 + II (:^2V3^II 2,-5/2 + l|3^ll2,o). (3.39) 
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Now the separate terms can be easily bound 



II C/liVll 2,-5/2 < l|C/lll2,-5/2l|J^lloo.O 
<C||3^||l^4,0 

<C(||y||4,0+||VJ^||4,-l) 

<c(iiy/y/^iu,o + iivy/*v:y^/*ii4,-i) 

< C(||3;||^fo^||3;||^/oV llVJ^II^f^illVJ^II^f^i) 



<cii3^ii;/2'oii3^ii2:'2,o 

<e||3^ll2,2,0+3ll3^lll,2.0, 

e 



where the last line conies from Young's inequality. 
Similarly, we have 

l|t/2V3^ll2,-5/2 <c||C/2||6,-3/2l|V3^ll3,-l 
<c||C/2||l,2.-3/2l|V:y||3,-l 

< C||V(3^)^/'v(3^)'/'||3,-i 

< C||V3^||^/!i||V3^||^/!i 

<c\\ijy\\y,\\\ijy\\l'l_, 

<e||3^ll2,2.o + -5113^111,2.0, 
e 

where we have switched from c to C to indicate the constant's dependence on || (5, A, tt, £)||jr. 
Combining these estimates with p.39p and applying (|2.7p to ||3^||i 2 '^^ have 

113^112,2,-1/2 < c(||^^(3^)||2,-5/2 + I|3^ll2,0)- (3-40) 

By construction, the adjoint operator F* has the same structure and thus also satisfies an estimate 
of the form (13.401) . In particular, this implies ker(F*) C W'^y^- Take a sequence G kcr(i^*) 

with ||A'„||2,2,-i/2 < 1 and we have WX^ - -^m||2,2,-i/2 < C'||A'„ - X^W^.q and by passing to a 
subsequence and applying a weighted version of the Rellich compactness theorem (see [lOj , Lemma 
2.1), the closed unit ball in ker(i^*) is compact and thus ker(i^*) is finite dimensional. The same 
reasoning tells us that ker(F) is also finite dimensional, thus there is a closed subspace Z such that 
W"^"^!^ = Z (B ker(F). To show F has closed range, it will suffice to prove 

113^112,2,-1/2 < C||F(J^)||2,_5/2, for all yeZ. (3.41) 

If (|3.4ip were not true, we could take a sequence 3^„ in Z with 113^,1112.2,-5/2 = 1 and ||-F(3^„)||2.-5/2 ^ 
0. Then we can again pass to a subsequence converging in Lq and (|3.40p implies converges to 
some y in Z O ker(i^), which would be a contradiction. Now we have ker(F) © coker(F) = Af* . 
Clearly Ran(i^) C Ran(£'$), so at most Ran(Z34>) differs from Af* by a finite dimensional closed 
subspace, and since Ran(_D$) ~ Af* , is surjective. 

□ 
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4 The Hamiltonian 



It is well known, that in order to generate the correct equations of motion, the first variation of the 
Hamiltonian density must be of the form 

5H = X ■ 5q + Y ■ 5p, (4.1) 

where q and p are the canonical position and momentum respectively. Hamilton's equations can 
then be read off as 

where t is the time parameter. 



In the framework of general relativity we need to make precise what we mean by 'time'. We 
interpret t as the flow parameter of the (yet to be specified) lapse-shift vector field on the spacetime. 
See [T3] for a detailed discussion on this. 

In the Einstein- Yang-Mills case, we have a vector field ^ on the bundle *P generating the evolu- 
tion, in place of the usual lapse-shift vector field. This corresponds to both a choice of coordinates 
on ''A^ and a choice of gauge. One may interpret the flow of ^ as simultaneously evolving the data 
through time while continuously changing the gauge. See [2J for details. 

In the Einstein- Yang-Mills case, the canonical position and momentum are {g,A) and (7r,e) 
respectively. In order to generate the correct equations of motion, we should expect to write the 
first variation of the Hamiltonian density as 

5H{g, A, TT, e) =X^ (g, A, n, e; ff'Sg,^ + X^{g, A, -k, e; (4.3) 

-f yi(5, A, ^, e- 0^J5n"' + Y^ig, A, tt, e; C) ■ fcl, (4.4) 

or equivalently 

DH^g^A,.,e) ■ {hApJ) = (^l,^2,i"l,>'2)l(s.A.,s) • {hAp,f)- (4.5) 

Notice that the Hamiltonian must have some dependence on the direction ^ in *P, in which we are 
to evolve the data. 

For the remainder of this paper, it will be convenient to write a point in the phase space as 
G — {g,A,-K,e) G F and a tangent vector Z — {h,b,p,f ) e TqJ^. With this notation, the usual 
ADM Hamiltonian (with Yang- Mills fields) is given by 

n^''^\G;0 = - I C'^ciG), (4.6) 

JM 

this is the pure constraint form of the Hamiltonian. 

It will be shown, if the boundary terms which were cast out when defining D<^* do indeed vanish, 
then we have 

DU^^'\0-Z^ - f D<S^UO-Z, (4.7) 

JM 
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and therefore the Hamihonian density is of the form (|4.5I) . We will see however, that if £, is not 
asymptotically zero, then this will not be the case; we will discuss this point later. 
Hamilton's equations now become 



d_ 

dt 



where 



J 






-1 




1 

1 



-10 



(4.8) 



(4.9) 



is the natural symplectic structure on J^. Equation (|4.8|) motivates Moncrief 's result [12] equating 
elements of ker(I?$g) with spacetime Killing vector fields. We thus refer to elements of ker(iD$Q) 
as generalised Killing vectors and if ^ corresponds to a stationary Killing field, then we call G = 
(g, A, TT, e) generalised stationary data. 



Proposition 4.1. The map H 



ADM 



FxAf - 



is a smooth map. 



Proof. The smoothness in G follows from the smoothness of $. We have jH'^^*^ (G; ^) | 



$||2._5/2, that is, 11^^^^ is bounded and linear in ^. 



l/2l 



We next establish the validity of equation ()4.7p . 



le$lll,-3< 
□ 



Theorem 4.2. For all i G W_( 



II' 



DUi^''\i) ■ Z = - [ D^UO-Z, 

J M 



for all Z e TqJ- . 

Proof. This is equivalent to the statement that the formal adjoint of D^q, given in section [3l 
is indeed the adjoint. We simply must demonstrate the boundary terms at infinity arising from 
integration by parts do indeed vanish. These boundary terms are given by 



Z-D^UO-eD<fGc.{Z) = 

V'^((7V(V,tr<,/i - V-'/i,,) + V\N)h,^ 



tTghV,iN))^-2X^p,^+Vy,, 



(4.10) 



The boundary terms have been expressed as two separate divergences corresponding to their decay 
rates at infinity - this distinction will be important later. Note, these divergences do indeed make 
sense as boundary integrals at inifinity in the usual (trace) sense, (see [5] - Lemma 4.3 and Lemma 
4.4). Formally, we consider an exhaustion of M by compact sets A^^, with smooth boundary, and 
take the limit of the boundary integrals as fc — )■ cxd. For convenience, we choose the exhaustion to 
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be euclidean balls near infinity and consider the limit of integrals on spheres. 
Lemma 4.4 of [5] gives us the estimate 

/ k|<c/R||7.||i,2,-3/2:A,, (4.11) 

where Sji is the sphere of radius R centred at zero, and Afj is the region bound between Sf; and 
52^. For simplicity, let us denote by V''BI and V';B,f , the first and second divergences in (|4.10|) 
respectively. is a collection of terms of the form q;/3^; where a G W'^y^, (i € ^"^^ 

^ G W^__\/2'' Note, g and g are bound, so we needn't consider the raising or lowering of indices 
in our estimates. 
Applying (|4.1ip . we have 

l^^l <c||a|L^sJ|e||oo:sJI/3|ll:S, 

Sr 

<o(l)||elloo:5jl/3||l,2.-3/2, (4.12) 

where we have made use of the fact a € W'^'^^^IM) C C''\M). In the limit as R tends to infinity, 

this integral vanishes and therefore B contributes no boundary terms. Note that this still holds if 
^ is only C° and bound. 

B"^ can be expressed as a collection of terms of the form ^t, where ^ £ ^-1/2 ^^^'^ ^ ^-3/2- 

/ \B'\<cmoo:sJrh.s^ (4.13) 
Jsr 

<o(l)||r||i,2,-3/2 (4.14) 
For the same reasons as above, the remaining boundary terms also vanish and therefore we have 

□ 

The necessity that ^ — >■ at infinity is twofold; not only does it ensure we have control on B^, 
it is also required for y^^^^^ to be well defined on J^, as <&(G') is not integrable for generic initial 
data. When ^ is taken to be asymptotic to some non-zero constant vector, it will be shown that 
the non- vanishing boundary terms, B^ , correspond to the first variation of energy- momentum and 
charge. This leads us to modify our Hamiltonian a la Regge and Teitelboim [20] . Before discussing 
this, we should make precise what we mean by "asymptotic to a constant vector" . 

Fix some € R'^'^ © g, which on some exterior region i?^, we may identify with e 
C°°(A° X TM X (g) a"), such that = 0. We now represent as G {A° x TM x g (g) A°), 
with = on E2B. and ^ = on Ej^. Obviously is not unique, however the difference 
between any two choices of is in Co" (A" X TTW X (0 (g A'')) c JV. This means that the space 

eoc+AA:=U:^-eooeAA} (4.15) 

is well defined. 
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Let us briefly digress to discuss the ADM energy-momentum and Yang- Mills electric charge. 
The ADM energy-momentum covector, ¥{g,T:) = {E,p.i) is usually given by 

lenE = / - d^gu)dS' (4.16) 

Wirp, = 2<j> TT.jdS^, (4.17) 

where the indices refer to some rectangular coordinate system at infinity. We use E to indicate the 
energy, to avoid confusion with the electric field. 

We also introduce the standard definition of Yang-Mills electric charge 

16^Q, 4 / E,,dS' = - / e,,dS\ (4.18) 
which clearly agrees with (|2.3p and the usual Maxwell total electric charge. 

It will be more convenient to work with integrals over A4 of divergences, rather than surface integrals 
at infinity. For a fixed G M.^'^ g, we define P in terms of its pairing with 

16<^Po(g) - / (el(V*V^5.j - Atr^g) + V^fl(V^g,, - V.tr^g)) ^ (4.19) 

167rCP.('r) - 2 / [I'Jj.i^l + V,L,) (4.20) 

JM ^ ' 

16^C^Pje) = 4 / [CJJ.K + K^^C^ , (4.21) 

JM ^ ' 

where indices are raised and lowered using the background metric, g. We know from |5j that 
P^ = (Po,Pj) is C°° on the constraint submanifold (if the source is integrable) and the definition 
is independent of g. Since i? is a density, the definition of P^j is clearly independent of °g and it is 
straightforward to check it is smooth. 

Proposition 4.3. // cr e L {tC \M) g) and {s,S,,a) € A^*, then ¥a{e) is a smooth function on 
the constraint submanifold, 

P„eC°°(C(s,5,a)). (4.22) 

Proof Fix some and choose a representative G Coo+-^- Since V(^;^) is compactly supported, 
we have 

^^\eM<\\L<7E\\, + \\Ei/u\i 

< \\L\Uma\\i + \\[A,E]\\,) + \\E\\,,.y2 

< llecollccXIklll + P(lll,2,-3/2l|i?lll,2,-3/2) + 11^112,-3/2- 

P„ depends linearly on E and is bounded, completing the proof. 

□ 

For ^ G + -^i we define the modified Hamiltonian, 

H«^(G;0 = 167rrF„- / C^c, (4.23) 

JM 
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recalling G — {g, A, tt, e). 

On shell, this new Hamiltonian now gives a value for some kind of total energy of the system. One 
should note that critical points of this Hamiltonian correspond to constrained critical points of the 
energy, with ^ acting as the (infinite dimensional) Lagrange multiplier. Unfortunately, neither of 
the terms in (|4.23p are well defined on all of however, it will be shown that the dominant terms 
cancel out. To see this, we define the regularised Hamiltonian 

n{G;0= f {L~-0-^+ f ^l(V'V^g-A(tr^g)v^-$o) (4.24) 
Jm j m 

+ / V^Cl ( V^5y - V, tr^ 9)^9+1 too (2V,7r^' - $ J 
Jm Jm 

+ f 27r'^ V,|^^. + / e^(4V,i?:-$J+4 / E^W.too- 
Jm Jm Jm 

We have combined the terms in (|4.23p and then separated them out into 7 integrals, each of which 
can be shown to be finite. 

Theorem 4.4. The regularised Hamiltonian defined by ^ . 24^ is a smooth functional on J- x (^^^ + 
Af). Furthermore, i/^ G + W'^f^j^, then for all G ^ J- and Z G TqJ- , we have 

DU^G.oiZ) ^ - I Z-D^U^). (4.25) 
Jm 

Equivalently, the regularised Hamiltonian generates the correct equations of motion. 

Proof. First we must establish boundedness and then smoothness follows from the same argument 
as Corollarv 13.21 

The first integral is easily bounded by ||^ — Cooll2,-i/2ll*5'll2,-5/2- The second and fourth integrals 
are bounded by Proposition 4.2 of 0- The 3rd, 5th and 7th integrals are bounded because V^qo 
has compact support, leaving only the 6th term; for this we note 

4V-i?,-$, =4[A,i?]„ (4.26) 

which is easily taken care of. This establishes the smoothness of we now prove the validity of 

Theorem 14.21 allows us to rewrite the variation of the first term in (|4.24p as Z ■ D^q{£^^ — ^); 
we will consider variations of the remaining terms separately. Consider the variation of the second 
and third terms, 

/ {<7\i'L{^' - V,tr^ h))^- Viili^ni,, ~ V, tr, h))^ (4.27) 
Jm 

+ V'iKjV'i"^ - trg hVdl)V9 - (h, b,p, f) ■ 
where the middle two terms in the above expression arise from the difference 

z-D^Uil)-il-D^Giz). 
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The third term in expression (j4.27p vanishes as V(^^) has compact support. The dominant terms 
in the first two divergences cancel, leaving us with a boundary term of the form ft,, ignoring 
factors of g. Note G leaving us a boundary term f /i, which is of the exact form of 
considered above and thus contributes nothing. Similarly, the variation of the fourth and fifth 
terms in (|4.24p give 

2/ |v,(f5,fep'=)-V,(elp})-Z-i?<l'*(G)(el)|. (4.28) 

The first and second terms in (j4.28l) give the boundary term [g — g)p, which is again of the form of 
above. Finally, the variation of the 6th and 7th terms in (|4.24l) give 

/ |-4V,(el/:) + 4V,(el/:) - Z ■ DKioitoo)} , (4.29) 

where the first and second terms here exactly cancel. Putting all of this together completes the 
proof 

□ 

5 The First Law 

It is well known that there is a strong analogy between the laws of thermodynamics and those of 
black holes. The first law is usually expressed by the in terms of differentials as, 

dm= —dA + fldJ + VdQ, (5.1) 
Stt 

valid for perturbations of stationary solutions. Here, the quantities m, k, A, i7, J, V and Q 
correspond to the mass, surface gravity, horizon area, angular velocity, angular momentum, electric 
potential and electric charge of the black hole, respectively. With our conditions for M and g, there 
will be no black hole present so this expression will reduce significantly. Also, interpreting V as 
the potential difference between the horizon and infinity, we expect to replace V with — l^o in this 
expression, leading us to 

dm + V^dQ = 0. (5.2) 

Theorem 15.21 and the subsequent corollary provide a proof of this similar to that of Sudarsky and 
Wald, as well as a converse statement conjectured in [22], which could not be rigorously shown 
without the Banach manifold structure given in Section [3] We prove that a solution satisfying this 
version of the first law must be stationary. Sudarsky and Wald also conjectured that a similar result 
should hold when an interior boundary is present, however it is likely that the precise form of (|5.ip 
be modified or the phase space include boundary conditions. A potential candidate for suitable 
boundary conditions are the conditions of an isolated horizon, as Ashtekar, Fairhurst, Krishnan 
and Beetle (see [31 |B] and references therein) have established a local version of the first law for iso- 
lated horizons. However this is beyond the scope of this paper and will be considered in future work. 

To begin, let us first quote the following generalisation of the method of Lagrange multipliers 
to Banach manifolds (Theorem 6.3 of |5]). 
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Theorem 5.1. Suppose K : Bi ^ B2 is a C map between Banach manifolds, such that DK^ : 
Ty^Bi — 7V(u)^2 is surjective, with closed kernel and close complementary subspace for all u G 
K ^(0). Let f G C^{Bi) and fix u ^ K ^(0); the following statements are equivalent: 

(i) For all v G kci DK^, we have 

Dfuiv) = (5.3) 

(a) There is X £ B2 such that for all w G 

Df^{v)^{X,DK^{v)), (5.4) 
where ( , ) refers to the natural dual pairing. 



We can now apply this to prove the main result. 

Theorem 5.2. Take G = {g, A, tt, e) G J" such that $(G) = {y, F,, 7a) G (a^ x T*M (g) x g* (g) 

-^6^ (fooj Coo) £ R^'^ (B g be fixed and define the energy functional E G C°°(C(s, 5, a)) by 

E{G)^Co^^{G), (5.5) 
The following statements are equivalent: 
(i) For all Z = {h,b,p, f) £ TqC{s, S,a) 

DEg{Z)=0 (5.6) 
(a) There is ^ G + W'^y^ ii'T ^ fl)) satisfying 

D^*G^ = 0. (5.7) 

Proof First we show (i) (m). For any fixed | G + ^-1/2' define /(G') = 'H{G';Cj for all 
G G J" and A' = $ - (s, S*, cr). With w = G, we now have the hypotheses of Theorem 15.11 Note 
also, TqC{s, 5, a) — ]iev{DKQ) and on TqC{s, S, cr), we have Df^ = WttDEq. So directly applying 
Theorem 15.11 there exists A G A/" such that for all Z G TqJ^, 

Dfg= I \-D^h{Z). (5.8) 

JM 

Inserting the definition of /, we have 

D<^1;{\) ^ ~D<^>*aii) (weakly). (5.9) 
Applying Theorem 13.31 we have 



D'S>UO = 0, (5.10) 



where C = X + ^e^^+ W_'^^^. We now have (i) (ii). 
To show (a) ^ (i), we simply must recaU that IGTrDEoiZ) = DndZ) for aU Z G TGC{s,S,cr) 
then from Theorem |H] we have mnDEaiZ) = DUaiZ) = -Jj^Z ■ D<^>a{^) = 

□ 
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Recall, a solution {N'^,V°') = ^ of D^q{£^) = corresponds to a generalised Killing vector N'^ 
and 'electric potential' V°' , representing an infinitesimal symmetry in the bundle. That is, evolution 
along integral curves of ^ in the bundle leaves the data fixed (cf. [2]). It was shown by Beig and 
Chrusciel[3, that if a Killing vector is timelike at infinity then it is asymptotically proportional to 

= rj^'^P^, where 77 is the Minkowski metric. It was further shown that, provided Tqq > T^qT^^ 
and P 7^ 0, C'^IP/i > for all future timelike vectors, C,. We will say the covector P^ is future timelike 
if ri'^'^Vi, is timelike and C^P^ > for all future timelike vectors, C,. 

Corollary 5.3. Suppose G G J-", ^(G) — (s, S", a) G and P^ is future timelike, then the following 
statements are equivalent 

(i) For all Z e TaC{s,S,a); 

Dma{Z) + V^-DQa{Z)=0, (5.11) 

where m = y—F^^ is the total (or "rest") ADM mass, V^^ £ g is the (fixed) Yang-Mills 
electric potential at infinity, and = E^^ds^ is the Yang-Mills electric charge. 

(ii) G is a generalised stationary initial data set with infinitesimal symmetry generator {N'^ , V°') — 
^, in the sense D^q{£^) = 0, and Ni!^ is proportional to . 

Proof. First we show (z) {ii). Choose = ——rj^'^F^, a future pointing unit timelike vector, 
such that N^P^ = m. For = (iV^,F^), we have DE = CDP^ = Dm + V ■ DQ and thus 
Theorem 15.21 gives us (ii). 

Conversely, (ii) implies the condition (ii) of Theorem 15.21 so we have DEq{Z; N^, V^) = for all 
Z e TqC{s, S,a). Since we have N^^ proportional to P*", we can rescale ^ such that N^F^ = m 
again, which completes the proof. 

□ 
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